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Using the recent thermodynamical study of isolated horizons by Ghosh and Perez, we provide a 
statistical mechanical analysis of isolated horizons near equilibrium in the grand canonical ensemble. 
By matching the description of the dynamical phase in terms of weakly dynamical horizons with 
this local statistical framework, we introduce a notion of temperature in terms of the local surface 
gravity. This provides further support to the recovering of the semiclassical area law just by means 
of thermodynamical considerations. Moreover, it allows us to study the radiation process generated 
by the LQG dynamics near the horizon, providing a quantum gravity description of the horizon 
evaporation. For large black holes, the spectrum we derive presents a discrete structure which 
could be potentially observable and might be preserved even after the inclusion of all the relevant 
transition lines. 



Introduction. In order to fully understand and ex- 
plain the semiclassical result of Hawking's calculation 
regarding the emission of thermal radiation from a black 
hole, it is strongly believed that a theory of quantum 
gravity is necessary in order to take into account the 
quantum structure of the horizon geometry at the Planck 
scale. 

The issue of the dependence of Hawking thermal spec- 
trum derivation on the behavior of the fields at arbitrarily 
high frequencies has been analyzed, for instance, by Ja- 
cobson P] and Unruh [3], with the first one speculating 
that, due to the presence of a short distance cutoff, the 
spectrum may undergo modifications. 

The possibility to probe the Planck scale structure of 
black holes with observations at much bigger wavelength 
has been conjectured by Bekenstein and Mukhanov [4]. 
Developing earlier proposal of theirs regarding black 
holes as discrete quantum systems, the authors start from 
the assumption that the area of the black hole be quan- 
tized and, due to the relation between black hole surface 
area and its mass, they study the emission amplitude re- 
lated to the jump from one quantized value of the mass 
to a lower one, in analogy with atoms behavior. Due to 
the discreteness of the area eigenvalues, the full emission 
spectrum is given by a set of spectral lines at frequen- 
cies multiple of a fundamental one of order of the black 
hole surface gravity, whose envelope reproduces Hawk- 
ing's thermal spectrum. This would imply a drastic de- 
parture from Hawking's semiclassical result, as empha- 
sized by Smolin [5]. 

However, Rovelli et al. have shown [B] that, replac- 
ing the naive ansatz of j4] for the area spectrum with the 
one computed from LQG, the spacing between the energy 
levels decreases exponentially with the black hole mass 
and, henceforth, for a macroscopical black hole, the spec- 
tral lines are very dense in frequency, recovering in this 
way a continues spectrum (see, however, [7] for a more 
detailed analysis based on the LQG area spectrum and 
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[S] for a numerical investigation applied to microscopical 
black holes). 

The 'atomic' picture of a quantum black hole adopted 
in was further exploited by Krasnov ^ within the 
LQG framework, soon after the first exciting ideas [TU] 
concerning the derivation of black hole entropy from the 
d.o.f. of the horizon quantum geometry started to cir- 
culate. Those ideas led to the introduction of a local 
definition of black hole suitable to a quantum gravity 
treatment, from which a microscopic origin of black hole 
entropy was derived. More precisely, through the quasilo- 
cal notion of Isolated Horizons (IH) [llj , able to capture 
the main physical properties of black holes event hori- 
zons, the LQG techniques were applied to space-time 
regions containing such an horizon, leading to a model 
[12l|T3] where the quantum black hole d.o.f. are described 
by a Chern-Simons theory on a punctured two-sphere. 
The punctures represent the quantum excitations of the 
gravitational field, as described by the LQG kinematics. 
Counting the dimension of the Hilbert space of such a 
theory living on the horizon leads to an entropy in agree- 
ment with the Bekestein-Hawking result, once the large 
area limit is taken. 

Krasnov's original idea was to use this quantum me- 
chanical description of black hole states to study the 
emission and adsorption spectrum of quantum black 
holes. By analogy with the gas of atoms and the use of 
Fermi's golden rule, he analyzed the intensity of emission 
lines corresponding to a single puncture transition from 
spin j to j' . The spectrum he finds shows a thermal 
character, even though it presents a discrete structure. 
However, in his analysis no dynamical process responsi- 
ble for the puncture transition is taken into account. In 
particular, in the derivation of the spectrum he considers 
a single atom transition approximation; this leads to the 
awkward feature of removing the line 1/2 0, which, for 
macroscopic black holes, would represent the most likely 
transition. Moreover, the statistical mechanical frame- 
work of [31 lacks a clear connection with the usual energy 
canonical ensemble, which is the fundamental ingredient 
for the derivation of all the thermodynamic functions (see 
[T3] for a recent analysis of the area canonical ensemble) . 
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In this letter, we want to investigate further and more 
in detail the analogy between a quantum horizon with its 
punctures and a gas of particles by introducing the main 
ingredients for a grand canonical ensemble analysis[^ The 
basic idea is to regard the bulk and the horizon as form- 
ing together an isolated system. The two subsystems are 
considered separately in thermal equilibrium; then, at 
some point, a weakly dynamical phase takes place and 
they interact with each other. This local interaction al- 
lows for the possibility of exchange of energy and particles 
between the two. After such a change of thermodynamic 
state has taken place, the two subsystems go back to a sit- 
uation of thermal equilibrium. This picture will be made 
more precise in the following, where we will concentrate 
only on the spherically symmetric case. However, let us 
at this point clarify the framework we are working in: no 
background structure is introduced at any point, we will 
work in the quantum gravity regime; no matter is go- 
ing to be coupled to gravity; the radiation spectrum we 
will derive is related entirely to emission of quanta of the 
gravitational field due to dynamical processes described 
by the LQG approach. 

Now, in order to provide the tools for a thermody- 
namical study of the system described above, one has to 
introduce a notion of local energy for the horizon. As 
it is well known, in the context of IH there is no unique 
notion of quasilocal energy function. This is due to the 
fact that there can be radiation in space-time outside the 
horizon and, therefore, no unique time evolution vector 
field can be single oul|j Nevertheless, in [121 EZj a local 
first law for isolated horizons has been recently derived, 
whose uniqueness can be proven once a local physical 
input is introduced. Interestingly, the notion of energy 
associated to the isolated horizon for a preferred family 
of local observers ^, at a proper distance (. from the hori- 
zon, turns out to be proportional to its area A, namely 



E = 



SttG 



A, 



(1) 



where R— 1/1 + o{£) represents the local surface gravity 
measured by the locally non-rotating stationary observer 
[T7] . This result sets the basis for the study of the 
quantum statistical mechanics of IH started in [16] . 

Entropy. Let us now first concentrate on the deriva- 
tion of the entropy of the gas of punctures (see [18] for the 



original microcanonical derivation and [TH] for a recent 
review). By working in the grand canonical ensemble — 
which represents the physically most suitable framework 
to describe the horizon+bulk system — , it can be shown 
how the Bekenstein-Hawking semiclassical entropy can 
be recovered only through thermodynamical considera- 
tions. Moreover, the description of the radiation process 
in the second part of the paper justifies the interpretation 
of the local notion of surface gravity introduced above as 
a temperature, which is a fundamental ingredient to re- 
cover the semiclassical area law (see below). In this sense, 
the result of the second part of the paper puts on more 
solid ground the recent derivation of [T^]. This section 
simply presents a more detailed derivation of the IH en- 
tropy in the grand canonical ensemble already performed 
in [16 . The original part of the paper is contained in the 
next section. 

The grand canonical partition function for the gas of 
punctures is given by 



^(/3)= ^z^Z(/3,iV)^ 



(2) 



N=0 



where z = exp{/3fi), fi being the chemical potential, N 
the number of punctures and Z(l3, N) the canonical par- 
tition function given 



Zip,N) 



1)^ 



(3) 



{Sj} j 



In the previous expression, {sj} represents a given quan- 
tum configuration where Sj punctures carry the spin-j, 
for all possible value of j, and such that J^j = -^i while 
Ej is given by the scaled IH area spectrum through ([I]), 
namely 



H\{s,})^y2s,E,\{s,})- 



j 



G 



sjV7(7TT)|{sj}). 

(4) 

Notice that, so far, (5 can be regarded simply as the inten- 
sive parameter conjugate to the energy ([T]) and its form is 
left unspecified. By means of the multinomial theorem, 
the grand canonical partition function ^ takes the form 



N 



oo / 



N=Q 



^ See |15| for a previous grand canonical ensemble treatment of a 
gravitational system with a spatial boundary, where the Hamilto- 
nian used in the definition of the partition function is the quasilo- 
cal energy associated with the boundary. 

^ Note that the same ambiguity is present also in the case of Killing 
horizons, since the freedom to rescale the Killing field normal to 
the horizon is still present. However, in the case of asymptotically 
flat space-times admitting global Killing fields, the normalization 
of the Killing field can be specified in terms of its behavior at 
infinity. 



(5) 



Given the partition function, all the thermodynamical 
functions of the system can be computed from it. In 



^ Here we are assuming the punctures statistics to be bosonic. We 
will comment more on this at the end of the paper. 
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particular, for the entropy we get 
d ,1 



S = -p-' — i-logT) = pE- fi(3N + \og^, (6) 



where 



E = ~—\og3f and N = z— log ^ (7) 
op oz 

are respectively the mean energy and the mean number 
of punctures of the gas. If we compute the expectation 
value for the occupation numbers sj, namely 

- ' ^og^^^-^m^^^^ (8) 



it is immediate to see that equations ^ imply E = 
y^^- SjEj and N = Sj. Therefore, the grand canonical 
entropy ([6]) can be written 



as 



I3R 
8^ 



A - nl3N + log jT, 



(9) 



where the first term in the r.h.s is the microcanonical en- 
tropy computed from the dominant configuration of the 
occupation numbers {s^} [16j. The expression (j9| can be 
seen as the horizon entropy at a generic temperature /3 
for the preferred local observer hovering outside the 
horizon at proper distance as resulting from the (phys- 
ical) local version of the IH first law dE — R/ {?>TiG)dA, 
which implies the local notion of energy ([!]) il7ij. 

In order to recover the Bekenstein-Hawking result from 
the grand canonical entropy ([9]), we need to look at its 
form in the semiclassical regime. First, we take the large 
N limit; secondly, we need to introduce a physical input 
proper of the black hole geometry outside, but near, the 
IH. 

When the number of punctures is big, the second of 
gives 



l^N/{N+l) = z^{2]^l)e- 



from which 



-l//31og(^(2j + l)e-''^0- 



(10) 



(11) 



Notice that, by means of the expression ^ for the en- 
ergy levels, the sign and the vanishing of the chemical 
potential above become a function of 7; this may have 
important physical implications for different sectors of 
the value of 7, which, however, won't be investigated 
here. As pointed out in [T^, the chemical potential does 
not need to vanish in order to recover the semiclassical 
result, due to the quantum nature of the hair correction 
iV/i. For the gas of photons description of black-body 
radiation, the vanishing of the chemical potential is con- 
sistent with the requirement that the free energy of the 



system be a minimum, since the average total number 
of particles in the gas is indefinite. In our case, N is 
not fixed a priori either, but, due to the dynamical pro- 
cesses taking place near the horizon (as described below), 
a change in the number of puncture implies a change of 
the horizon energy and, therefore, /i can be different from 
zero. However, semiclassical consistency now implies a 
quantum modification of the first law, where a chemical 
potential term has to be introduced due to this hair term 
in the entropy. 

The last missing ingredient, in order to recover 
Bekenstein-Hawking area law [TJ |20] , is now the identifi- 
cation of the radiating temperature /3 with the inverse of 
the horizon surface gravity. Let us recall that this identi- 
fication made sense only after Hawking's surprising result 
that the black hole could radiate; before that, it only had 
the status of a resemblance or a coincidence suggested by 
the equations form, since the black hole temperature was 
expected to be zero. 

In the second part of the paper, we will argue that 
quantum geometrodynamical considerations could gen- 
erate a radiation process from the quantum horizon, for 
which the local surface gravity k can be properly inter- 
preted as a temperature. Despite being of a completely 
different nature (respect to the Hawking effect) and not 
specifying the proportionality coefficient, this process al- 
lows us to introduce a notion of temperature for the quan- 
tum IH. One can now fix this relationship between tem- 
perature and surface gravity using the semi-local (and 
semi-classical) argument of J16] . Namely, if we assume a 
stationary near-horizon geometry and we use the Unruh 
temperature /3k = 27r/fi, for our local accelerated observer 
^, the entropy expression ^ gives exactly S = A/{A£p), 
at the leading order. 

Radiation. We now switch our attention to the ra- 
diation process. The topology of the null surface A rep- 
resenting the IH is assumed to be x K. This means 
that, at 'a given time', the horizon surface is given by the 
two-sphere at the intersection between A and a Cauchy 
surface representing space. The quantum space geometry 
is described by the polymer-like excitations of the gravi- 
tational field encoded in the spin networks states, which 
span the kinematical Hilbert space of LQG. Some edges 
of those states can now pierce through the horizon sur- 
face, providing local quantum d.o.f. accounting for the IH 
entropy, as shown above. Dynamics in the bulk is imple- 
mented by the Hamiltonian constraint H[N] = 0. This 
reflects the fact that, in a diffeos invariant theory, the 
canonical Hamiltonian generates evolution in the param- 
eter of the action which is unphysical, i.e. the evolution 
is pure gauge. This freedom in choosing the (unphys- 
ical) evolution parameter is reflected in the possibility 
to rescale the time evolution vector field by the lapse 
N. Then, in the quantum theory, the projection opera- 
tor into the kernel of the Hamiltonian constraint is con- 
structed by integrating over all the possible value of the 
lapse N ([Sum]). 
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Despite the lack of a completely well defined quanti- 
zation prescription for the Hamiltonian operator, some 
general properties are satisfied by all the several propos- 
als present in the literature. Namely, H acts locally at 
the vertices of the spin networks states and its action 
changes the spin associated to some of the edges attached 
to the given vertex. In order to be more quantitative, in 
the following we will consider Thiemann's proposal |23| . 
Therefore, if we concentrate on the Euclidean constraint, 
for simplicity, its action on a 3-valent node having two 
edges piercing the horizon can be graphically represented 
as 





where the holonomies entering the regularization of H 
are taken in the fundamental representation (the funda- 
mental vertex on the right provides the covariant rep- 
resentation of the Hamiltonian operator action). The 
action depicted above can be interpreted as an absorp- 
tion/emission process of quanta of the gravitational field 
by the quantum horizon. 

In order to study the radiation process, we need to 
'perturbate' the equilibrium states represented by the IH 
configurations by turning on some weakly dynamical ef- 
fects near the horizon, till another equilibrium configura- 
tion state is reached again. In this way, a physical process 
can causes the transition from one equilibrium state to a 
nearby one. The horizon goes through a non-equilibrium 
phase resulting in a small change of its state. For a large 
black hole, this process is expected to be quite slow. 

In order to describe this evolution process, we need first 
to deparametrize the system near the horizon, i.e. we 
need to locally single out one of the partial observables 
to play the role of time. In this way, we can use the 
Hamiltonian operator in the bulk near the horizon to 
define evolution of the boundary quantum states with 
respect to this observable. 





FIG. 1. After the collapse of a spherical star, the IH Ai is 
created. Later, in the dynamical phase H, flux of gravita- 
tional and matter energy goes through the horizon, before a 
new equilibrium state A2 is reached again. 

Let us now show how the contact between the Dy- 
namical Horizons framework developed in |24l 125) and 



the thermodynamical description of [ini HZ] allow us to 
single out such a time evolution vector field. Within a 
framework advocated by Ashtekar, the idea is to describe 
the evaporation process by interpreting this intermediate 
phase of horizon state evolution as an extension of the 
quantum horizon geometry from isolated to dynamical 
horizons, where a proper quantum notion of gravitational 
energy — induced by the quantum theory of geometry in 
the bulk — replaces the classical one introduced in [251 f2^ . 

In |25] , a dynamical process version of the first law was 
derived from an area balance law relating the change in 
the area of the dynamical horizon to the flux of matter 
and gravitational energy. In the absence of matter, this 
could be written as 



/ r2 _ n_\ _ 
V2G 2g) ~ 



1 

167rG 



(13) 



H 



where H is the portion of dynamical horizon bounded by 
the 2-sphere leaves 5*1 , S2 of radius ri , r2 and Hg is the 
bulk term in the Hamiltonian, which can be written in 
terms of pure geometrical quantitiet]^ The infinitesimal 
version of ( 13 ) gives a first law for dynamical horizons, 
namely 



SttG 



dA = dEr 



(14) 



where = (2r)^^ is an effective surface gravity. 

First of all, let us notice that, given the spherical 
symmetry of the horizons studied here, the classical no- 
tion of gravitational energy defined by the r.h.s. of ( 13 ) 



would give a vanishing flux across the dynamical horizon. 
Of course, this doesn't mean that also at the quantum 
level there cannot be a flux of energy across the horizon. 
Hawking radiation taught us exactly the opposite. 

In |25| it is shown that, for a dynamical horizon, the 
first law (14) can be generalized to Rf / {%nG)dA — dEf, 
where Hf — (df /dr)Kr, with /(r) an arbitrary function of 
the radius r encoding the freedom to rescale the vector 
field — with respect to which the notion of energy is 
associatecPj— , i.e. = A^/x" = (df /dr)S,r, where Nf is 
a permissible lapse associated with the radial function /. 
The function /(r) encodes the dynamical nature of the 
previous version of the first law, since r plays the role of 
time along the dynamical horizon; but it also represents 
an ambiguity in the description of the horizon dynamics. 

The passage from a weakly dynamical horizon H to 
a weakly isolated horizon A is a well defined procedure 
described in |25, . The only ingredient missing in our con- 
struction is the matching of the surface gravity entering 
( 14 1 and the one appearing in the local notion of energy 



^ Notice that the I.h.s. of ( |13[ | is obtained from a boundary term 

in the Hamiltonian, as shown already in |27j . 
^ In the case of IH, this corresponds to the freedom to rescale the 

null normal vector field i and therefore the surface gravity, by 

a positive constant c. 



5 



Q. This last step can now be carried out by means of 
the rescaling freedom described above, using the local 
universal first law of [T7] to solve this ambiguity. 

For our preferred family of observers located right out- 
side the horizon, a natural choice of the field ^/ would 
be given by / — £, with £ being the proper distance from 
the horizor{^ If we use the fact that, in standard co- 
ordinates for the Schwarzschild metric, for an observer 
at r = 2M + e, the proper distance from the horizon is 
i = 2V2Me and set /(r) = £, we get an effective surface 
gravity matching exactly the local surface gravity enter- 
ing the energy expression (IT]), i.e. kj = ^/£V\ Therefore, 
the local first law for black nole thermodynamics derived 
in [ini \17\ is compatible with the framework of [25] if we 



single out, among the permissible time vector fields for 
the local observer ^, the one defined by = Nix"" ■ 

Moreover, if we interpret e as a parameter controlling 
the slow evolution of the horizon area, i.e. ^^rg ~ e, the 
previous expression for the surface gravity Rf matches 
the requirements of [25] for the definition of slowly evolv- 
ing horizons. In this case, the local surface gravity en- 
tering the relation ([T]) would be approximately constant. 
Henceforth, one can regard the horizon as making con- 
tinuous transitions from one equilibrium state to another 
(see FIG. [T]) and the geometrical surface gravity R can 
be interpreted as a (slowly varying) parameter control- 
ling this process. In this sense, it seems very natural 
to interpret R as the temperature of the radiation emit- 
ted during this transition phase. This strengthens the 
connection between thermodynamics and our statistical 
analysis and clarifies further the entropy derivation per- 
formed in [TO] . 

At the quantum level, the local notion of energy ([!]) 
acquires a definition in terms of the LQG area operator: 
this is the fundamental ingredient for the construction 
of the quantum statistical mechanical ensemble above. 
Since the Hamiltonian operator acts on spin network 
states and modifies them, we can now use the local de- 
parametrization of the system described above to inter- 
pret the bulk quantum dynamics near the horizon as a 
generator of boundary states evolution. In other words, 
we can now make sense of the r.h.s. of ( 14 ) at the quan- 
tum level by means of the LQG dynamics. The conse- 
quent jump to a different horizon area operator eigenstate 



® In the generalization of the area balance law introduced in 1251 . 
the function / is assumed to be constant in each 2-sphere cross 
section of the dynamical horizon H. In this way, the permissi- 
ble vector field is such that the flux of gravitational energy, 
associated to it, across a portion AH on the horizon relates the 
two constant values the function / assumes on the fixed cross- 
sections bounding AH. Since a change in the horizon radius im- 
plies a change in the proper distance of the local observer from 
one cross-section to another, our choice f = £ is compatible with 
the generalization of |25l . 

^ The same result can be obtained without referring to any global 
structure by using the local extension of the null normals x" to 
the vicinity of the horizon, as constructed in [TTj and used to 
define the preferred local observer 0". 



would thus result in an horizon energy variation. In this 
way, quantum dynamics would allow for an energy flux 
across the horizon, providing a quantum notion of weakly 
dynamical horizon. 

More precisely, the action of the Hamiltonian operator 
on nodes close to the horizon (see (12)) can be used, in 
the formalism of to define a physical scalar prod- 
uct to interpret as transition amplitudes between one 
boundary state to another, in the physical time evolution 
parameter defined by Ni. In general, given a quantum 
boundary state |{Sj}) on a given 2-sphere cross-section, 
the action of the Hamiltonian operator on the vertices 
near the horizon will evolve it in a state |{sj'}) through 
the change of spin of some of the punctures in the initial 
state The probability of the transition between 

an initial state |{sj}) and a final state |{sj'}) differing 
only at a given vertex v, in an infinitesimal interval of 
time, is given by Pjj' 



where 



^jj' — {{^j'}\{^j})phys 



= -im{{s,-}\H,\{s,})+o{N!). (15) 

Let us explain the last passage in the previous equation. 
In order to define a diff-invariant physical scalar prod- 
uct, the action of H is defined on states of J^aux, hence 
picking up a dependence on the regulator (i.e. the po- 
sition of the link added); however, this dependence fi- 
nally drops out in the integration over TV. More pre- 
cisely, one can insert inside {■\-)phys an integration over 
the diffeomorphisms and this confines the remaining diff- 
dependcnt quantities in the arguments of the functions 
N{x); the diff-invariant nature of the DN integral then 
eliminates any residual dependence. In the case of an 
internal boundary though, the situation is more subtle. 
In particular, IH boundary condition are such that only 
diffeos tangent to the horizon are to be regarded as gauge 
transformations in the Hamiltonian framework [TT] ; these 
diffeos are compatible with the 'gauge fixing' performed 
in the evaluation of (15). This gauge-fixing, which en- 



codes the local deparametrization of the system due to 
the choice of a preferred ('area') time variable described 
above, could be performed by inserting a delta function in 
the measure of the DN integral which singles out the con- 
stant value Ni for the scalar function N{x) at the vertices 
near the horizon. This procedure could be interpreted as 
the insertion of a local observable Oi — 5{N, Ni{v)) in 
the path integral representing our physical time-evolution 
variable. 

Let us now concentrate on the description of the radi- 
ation process. In defining Hjj> , we will use the matrix 
elements of the Hamiltonian operator obtained in [28], 
where only the Euclidean term, of Thiemann's ver- 
sion |23| of the constraint is considered. Eventual modi- 
fications of the radiation spectrum induced by the inclu- 
sion of the 'Lorentzian' term are left for future study. 

At this point, an important observation concerns the 
local nature of the dynamical process near the boundary 
in relation with the action of the Hamiltonian operator. 
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The terms generated by the action of H'^ on a generic 
3-vertex near the horizon are of two different kinds in 
our context: two of them create new punctures piercing 
through the horizon and one (that shown in (12)) does 
not. In the first case, the horizon area can only increase 
or at most remain constant (no radiation) . The inclusion 
of those terms creating new punctures though would spoil 
the local characterization of the interaction between the 
bulk and the horizon, by introducing an arbitrary de- 
pendence of the transition amplitudes on the spin of the 
bulk links adjacent to the interaction vertex. Moreover, 
since the new puncture created by these terms has to 
be co-planar to the two it is attached to (one piercing 
the horizon and one in the bulk), the entropy calcula- 
tion for the new equilibrium configuration may lead to 
divergences, due to the fact that one can no longer trace 
over the bulk d.o.f. without affecting the horizon state 
and the moding out of the diffeos tangent to the horizon 
becomes non trivial. 

For these reasons, we only allow the action of H'^ that 
doesn't create new punctures. Notice that this term cor- 
responds to the absorption inside the horizon of a quan- 
tum of the gravitational field (the new link created by 
ff^) associated to an emission process, due to the jump 
of the horizon state to a lower area eigenvalue. This 
dynamics is reminiscent of the heuristic picture of Hawk- 
ing radiation where one anti-particle is absorbed by the 
black hole while one particle escapes to infinity. In this 
sense, the original idea fP that Hawking had in mind, 
to describe the black hole evaporation, is realized in our 
context. 



Given these considerations, the Hjji in (15) corre- 
sponds to the transition amplitude where two of the 
punctures piercing the horizon have jumped to a differ- 
ent energy level. Since we want to study the emission 
process, we will now take into account only the jumps to 
a lower eigenvalue of the area (energy) operatoir] 

For the case of a macroscopic black hole considered 
here, only the first SU{2) Irreps are relevant for the 
horizon punctures, henceforth, we will consider only the 
cases where q = j and p = j,j + 1/2, j -I- 1, j -I- 3/2, i.e. 
r — 1/2, 1, 3/2 (r being the spin of the third link at the 
3-vertex acted upon not piercing the horizon). Higher 
values of r would not be relevant for the spectrum of a 
large black hole; this is a feature of the action of H when 
no new punctures are created, as observed above. We 
now choose the holonomies in to be in the funda- 
mental representation, implying that q can only jump to 
q-1/2 while p top ±1/2. 

Another source of ambiguity is represented by having 
two possible orderings of H'^, as defined by Thiemann. 
More precisely, the volume term inside H'^ can act before 



* This breaking of unitariety can be justified by the fact that we are 
looking only at a part of the system (the near-horizon dynamics), 
independently of what happens in the bulk. We thank C. Rovelli 
for pointing this out. 



or after the creation of the new link. Here we will choose 
the ordering in which the volume operator acts before the 
punctures have jumped in order to have a non vanishing 
probability for the transitions involving the jump 1/2 —> 
for one of the punctures. This will be our choice for 
the regularization of the Hamiltonian operator (we will 
comment below on the ambiguity of choosing different 
representations). 

Notice that, with this 'two atoms' approximation mo- 
tivated by the action of H'^, the jump j — >■ j — 1/2 in- 
volves always an even number of punctures; in this way, 
the gauge symmetry of the IH is preserved, as an equilib- 
rium state is reached again, and the unpleasant feature 
of removing the 1/2 — >■ transition from the spectrum is 
not present anymore. 



In the expansion ( 15 ) of the transition amplitude, we 
have kept only the first order term in Ni, corresponding 
to a single action of the Hamiltonian operator at v; this 
encodes at the quantum level the weak time dependence 
of the horizon in this intermediate phase between two 
equilibrium configurations. The contemporary action of 
Jf^ on more than one node would produce a bigger num- 
ber of lines in the spectrum, allowing for the possibility 
that the discrete structure of the spectrum derived below 
get smoothen, resulting in an effective continues spec- 
trum. However, there is the possibility that these lines 
would be significantly suppressed by the higher powers 
of the transition amplitude \N£{H^)\, preserving in this 
way the discrete structure. Further investigation of this 
feature is crucial for a possible observation of Planck's 
scale effects. 

The H'^ matrix elements A{2q, e; 2q, e; 2r) correspond- 
ing to the jumps 2p — >■ 2p + e, 2q — >• 2q + e, with e, e = ±1, 
are explicitly computed in [28]. 

We now have all the ingredients to derive the spectrum 
form, for our local observer ff, of the radiation emitted 
by the quantum horizon in the process — ?■ |{sj'}). 
Using the analogy with the gas of particles, the eigen- 
values (of the operator corresponding to ([!])) Ej repre- 
sent the energy levels of the system, which are occupied 
by Sj quanta of the gravitational field. In our picture 
then, the punctures play the role of the photons in the 
Bose-Einstein approach to the investigation of black-body 
radiation and Ej is the analog of a photon energy level 
hujp. According to the emission process described above, 
the energy of the quantum of radiation emitted after the 
action of H'^ on two punctures with spins p, q is given by 
A^^ = Ep — Epj. 1 + Eq — Eq_ 1 . Therefore, the intensity 
of the spectrum lines is given by 



Ijr = SpSg\NeA{2p, ±1; 2q, -1, r)\''A 



P9' 



(16) 



where, with the convention introduced above, q = j, j < 
p < j + r and j,r = 1/2, 1,3/2; the occupation numbers 
are given by eq. 

The total energy emitted in the process is given by / = 
X; Ijr- In FIG. [2] we plot the energy distribution Ij^. as 
a function of /3 times the energy of the quanta emitted. 
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FIG. 2. Emission lines witli tlieir intensities and the thermal 
envelope. The normalization factor A'^(iV + 1) of the intensity 
is left unspecified since we are interested only in the form of 
the spectrum. 

formed by four different sets of spectral lines. The first 
two groups of lines correspond to the case e = +1, i.e. 
when p p+ 1/2, while the last two to the case e — —1, 
i.e. p—>p - 1/2. 

The envelope of the emission lines reproduces a ther- 
mal behavior. However, since the frequency gap between 
the different sets of lines is of the order Auj de- 
parture from a continuos spectrum for macroscopic black 
holes could be potentially observable. 

Conclusions & comments. We have used the sta- 
tistical mechanical framework introduced in [Tni [TT] to 
study the properties of IH at and near the equilibrium in 
the grand canonical ensemble. This setting seems to be 
the most appropriate, given the possibility of the horizon 
to exchange energy with the bulk system, seen as a large 
reservoir, and the indeterminacy in the mean number 
of punctures. By matching the description of the inter- 
mediate dynamical phase in terms of weakly dynamical 
horizons [551 US] with the local statistical description of 
IH in [TBI HZ] J we could introduce a notion of tempera- 
ture regulating this exchange in terms of the local sur- 
face gravity Rg. Moreover, this allowed us to single out 
a physical time parameter in terms of which describing 
the boundary states evolution. In this way, our analysis 
provides further support to the main result of {16| , where 
the semiclassical area law for IH entropy is recovered just 
by means of thermodynamical considerations. 

The kinematical Hilbert space of LQG provides a beau- 
tiful description of the quantum horizon d.o.f.. By means 
of the regularization and quantization prescription of the 
theory for the Hamiltonian constraint, we managed to de- 
fine a quantum notion of gravitational energy flux across 



the horizon. This allowed us to study the radiation pro- 
cess generated by the Hamiltonian operator, providing a 
quantum gravity description of black hole evaporation. 
For large black holes, the discrete structure of the spec- 
trum could be potentially observable. 

We conclude with a few comments. In analogy with 
the criticism of [Sj to the Bekenstein-Mukhanov effect, 
one could ask if the inclusion of other transitions be- 
tween area eigenvalues would make the spectrum in FIG. 
[2] effectively continuos. According to the LQG dynamics, 
one way of adding more lines would be to go beyond the 
two atoms approximations and consider the simultaneous 
action of H on more than one vertex; as observed above, 
however, those lines could be significantly damped by the 
transition amplitude (151, leaving the discrete structure 
unalterecF^l 

Another possibility would be to use the ambiguity in 
the holonomies representation present in the definition 
of H. For large black holes, only the first few Irreps 
beyond the fundamental one would be relevant for the 
analysis. As it can be easily checked, due to the selection 
rules implemented by the Hamiltonian operator, the in- 
clusion of the first lowest Irreps produces only a few more 
sets of lines whose frequencies, except for a couple, are 
completely shifted from (and with the same pattern of) 
the ones given by the fundamental representation (with 
a gap again of order ^ k) or exactly overlap to them. 
Therefore, even though the shape of the spectrum might 
change by including more transition lines, the LQG dy- 
namics might preserve its discrete structure, providing 
the possibility to observe a departure from the semiclas- 
sical scenario. However, further investigation in this di- 
rection is necessary and our derivation of the emission 
lines plotted in FIG. [2] should be regarded as a first step 
towards a fully quantum dynamical derivation of the ra- 
diation spectrum. 

While our attention here has been limited to Thie- 
mann's proposal for the regularization of the Hamiltonian 
operator, the same analysis could be applied to different 
regularization schemes, such as [29 . The implications of 
alternative quantization of the Hamiltonian operator for 
the spectrum could eventually be used to discriminate 
between the different proposals and, hopefully, to reduce 
the number of ambiguities left in the definition of H. 

Moreover, the near-horizon dynamics could be imple- 
mented within the spin foam formalism. More precisely, 
given the increasing evidence for the spin foam amplitude 
to implement the Hamiltonian constraint in the canon- 
ical formalism 1301 . one could use the vertex amplitude 
to compute transitions probabilities between boundary 
states containing the right part of ( 12 ) as a subgraph. 



This could turn out to be a very useful application of 
the conjectured duality between the canonical and the 
covariant theory, with potential observational evidence. 



Here we have used the interpretation of the surface gravity as 
the radiation temperature discussed above and the fact that, in 

the large N hmit, fjR'flp/G ^ o(l). This would be the case, for instance, if N(N ~ o(l). 
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The analysis of these interesting issues is left for future 
study. 

Despite the implementation of the Hamiltonian con- 
straint to evolve the boundary quantum states, our anal- 
ysis is still an effective one, since based on an interplay 
between classical general relativity and quantum geome- 
try. It would be very interesting to apply the ideas devel- 
oped here to the framework of [5T] . where the quantum 
horizon d.o.f. are described within the full LQG formal- 
ism. This might provide further insights on the dynami- 
cal nature of the radiation process and could reveal unex- 
pected connections between three and four dimensional 
quantum gravity. Moreover, the analysis of [31' seems to 
indicate that the punctures must obey some kind of any- 
onic statistics; investigating the consequences that such 
a different statistics would have on the radiation descrip- 
tion could lead to important deviations from the bosonic 
assumption. 

Finally, on the relation with Hawking radiation. The 
two processes are quite different from each other. Hawk- 
ing effect is 'much less' dynamical |32| : it simply requires 
a curved background and a scalar field propagating on 
it: the gravitational d.o.f. are frozen; the emission pro- 
cess is intrinsically related to the definition of particle 
possible due to this switching off of gravity. In our anal- 



ysis, no background is ever picked up and no matter is 
coupled: gravity does all the job. It is the evolution 
of the quantum gravitational d.o.f. on the horizon, to- 
gether with its beautiful interplay with classical general 
relativity and thermodynamics, which is responsible for 
the energy emission. Moreover, due to this dynamical 
framework, the whole description of the process is local 
(no notion of future/past null- infinity is involved). 

On the other hand, as Hawking radiation is a genuinely 
quantum effect, the emission process described here is 
also purely quantum. As noted above, the classical first 
law ( 13 ) derived in [551 US] would give a vanishing fiux 
for spherically symmetric horizon; it's the evolution — 
described by the LQG dynamics — of fiuctuating quantum 
d.o.f. that generates radiation. 

The analysis presented here could be extended to the 
inclusion of matter, whose d.o.f. are also expected to live 
on spin networks edges [33^. Such an extension would 
provide a more appropriate set-up for a comparison with 
the Hawking effect. This is also left for future investiga- 
tion. 
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